BAR4R-PUB-06/073 
SLAC-PUB-12410 



Measurement of CP Violation Parameters with a Dalitz Plot Analysis of 

B± — D^-^K± 

B. Aubert, 1 M. Bona, 1 D. Boutigny, 1 Y. Karyotakis, 1 J. P. Lees, 1 V. Poireau, 1 X. Prudent, 1 V. Tisserand, 1 

A. Zghiche, 1 E. Grauges, 2 L. Lopez, 3 A. Palano, 3 J. C. Chen, 4 N. D. Qi, 4 G. Rong, 4 P. Wang, 4 Y. S. Zhu, 4 
G. Eigen, 5 I. Ofte, 5 B. Stugu, 5 G. S. Abrams, 6 M. Battaglia, 6 D. N. Brown, 6 J. Button-Shafer, 6 R. N. Cahn, 6 

Y. Groysman, 6 R. G. Jacobsen, 6 J. A. Kadyk, 6 L. T. Kerth, 6 Yu. G. Kolomensky, 6 G. Kukartsev, 6 D. Lopes Pegna, 6 
G. Lynch, 6 L. M. Mir, 6 T. J. Orimoto, 6 M. Pripstein, 6 N. A. Roe, 6 M. T. Ronan, 6 'B K. Tackmann, 6 W. A. Wenzel, 6 
P. del Amo Sanchez, 7 M. Barrett, 7 T. J. Harrison, 7 A. J. Hart, 7 C. M. Hawkes, 7 A. T. Watson, 7 T. Held, 8 H. Koch, 8 
B. Lewandowski, 8 M. Pelizaeus, 8 T. Schroeder, 8 M. Steinke, 8 J. T. Boyd, 9 J. P. Burke, 9 W. N. Cottingham, 9 

D. Walker, 9 D. J. Asgeirsson, 10 T. Cuhadar-Donszelmann, 10 B. G. Fulsom, 10 C. Hearty, 10 N. S. Knecht, 10 
T. S. Mattison, 10 J. A. McKenna, 10 A. Khan, 11 P. Kyberd, 11 M. Saleem, 11 D. J. Sherwood, 11 L. Teodorescu, 11 

V. E. Blinov, 12 A. D. Bukin, 12 V. P. Druzhinin, 12 V. B. Golubev, 12 A. P. Onuchin, 12 S. I. Serednyakov, 12 
Yu. I. Skovpen, 12 E. P. Solodov, 12 K. Yu Todyshev, 12 M. Bondioli, 13 M. Bruinsma, 13 M. Chao, 13 S. Curry, 13 
I. Eschrich, 13 D. Kirkby, 13 A. J. Lankford, 13 P. Lund, 13 M. Mandelkern, 13 E. C. Martin, 13 D. P. Stoker, 13 
S. Abachi, 14 C. Buchanan, 14 S. D. Foulkes, 15 J. W. Gary, 15 F. Liu, 15 O. Long, 15 B. C. Shen, 15 L. Zhang, 15 
E. J. Hill, 16 H. P. Paar, 16 S. Rahatlou, 16 V. Sharma, 16 J. W. Berryhill, 17 C. Campagnari, 17 A. Cunha, 17 

B. Dahmes, 17 T. M. Hong, 17 D. Kovalskyi, 17 J. D. Richman, 17 T. W. Beck, 18 A. M. Eisner, 18 C. J. Flacco, 18 
C. A. Heusch, 18 J. Kroseberg, 18 W. S. Lockman, 18 T. Schalk, 18 B. A. Schumm, 18 A. Seiden, 18 D. C. Williams, 18 

M. G. Wilson, 18 L. O. Winstrom, 18 E. Chen, 19 C. H. Cheng, 19 A. Dvoretskii, 19 F. Fang, 19 D. G. Hitlin, 19 
I. Narsky, 19 T. Piatenko, 19 F. C. Porter, 19 G. Mancinelli, 20 B. T. Meadows, 20 K. Mishra, 20 M. D. Sokoloff, 20 
F. Blanc, 21 P. C. Bloom, 21 S. Chen, 21 W. T. Ford, 21 J. F. Hirschauer, 21 A. Kreisel, 21 M. Nagel, 21 U. Nauenberg, 21 
A. Olivas, 21 J. G. Smith, 21 K. A. Ulmer, 21 S. R. Wagner, 21 J. Zhang, 21 A. Chen, 22 E. A. Eckhart, 22 A. Soffer, 22 
W. H. Toki, 22 R. J. Wilson, 22 F. Winklmeier, 22 Q. Zeng, 22 D. D. Altenburg, 23 E. Feltresi, 23 A. Hauke, 23 
H. Jasper, 23 J. Merkel, 23 A. Petzold, 23 B. Spaan, 23 K. Wacker, 23 T. Brandt, 24 V. Klose, 24 H. M. Lacker, 24 
W. F. Mader, 24 R. Nogowski, 24 J. Schubert, 24 K. R. Schubert, 24 R. Schwierz, 24 J. E. Sundermann, 24 A. Volk, 24 

D. Bernard, 25 G. R. Bonneaud, 25 E. Latour, 25 Ch. Thiebaux, 25 M. Verderi, 25 P. J. Clark, 26 W. Gradl, 26 

F. Muheim, 26 S. Playfer, 26 A. I. Robertson, 26 Y. Xie, 26 M. Andreotti, 27 D. Bettoni, 27 C. Bozzi, 27 R. Calabrese, 27 

G. Cibinetto, 27 E. Luppi, 27 M. Negrini, 27 A. Petrella, 27 L. Piemontese, 27 E. Prencipe, 27 F. Anulli, 28 
R. Baldini-Ferroli, 28 A. Calcaterra, 28 R. de Sangro, 28 G. Finocchiaro, 28 S. Pacetti, 28 P. Patteri, 28 I. M. Peruzzi, 28 '0 
M. Piccolo, 28 M. Rama, 28 A. Zallo, 28 A. Buzzo, 29 R. Contri, 29 M. Lo Vetere, 29 M. M. Macri, 29 M. R. Monge, 29 
S. Passaggio, 29 C. Patrignani, 29 E. Robutti, 29 A. Santroni, 29 S. Tosi, 29 K. S. Chaisanguanthum, 30 M. Morii, 30 
J. Wu, 30 R. S. Dubitzky, 31 J. Marks, 31 S. Schenk, 31 U. Uwer, 31 D. J. Bard, 32 P. D. Dauncey, 32 R. L. Flack, 32 
J. A. Nash, 32 M. B. Nikolich, 32 W. Panduro Vazquez, 32 P. K. Behera, 33 X. Chai, 33 M. J. Charles, 33 U. Mallik, 33 
N. T. Meyer, 33 V. Ziegler, 33 J. Cochran, 34 H. B. Crawley, 34 L. Dong, 34 V. Eyges, 34 W. T. Meyer, 34 S. Prell, 34 

E. I. Rosenberg, 34 A. E. Rubin, 34 A. V. Gritsan, 35 C. K. Lae, 35 A. G. Denig, 36 M. Fritsch, 36 G. Schott, 36 
N. Arnaud, 37 J. Bequilleux, 37 M. Davier, 37 G. Grosdidier, 37 A. Hocker, 37 V. Lepeltier, 37 F. Le Diberder, 37 

A. M. Lutz, 37 S. Pruvot, 37 S. Rodier, 37 P. Roudeau, 37 M. H. Schune, 37 J. Serrano, 37 V. Sordini, 37 A. Stocchi, 37 
W. F. Wang, 37 G. Wormser, 37 D. J. Lange, 38 D. M. Wright, 38 C. A. Chavez, 39 I. J. Forster, 39 J. R. Fry, 39 
E. Gabathulcr, 39 R. Garnet, 39 D. E. Hutchcroft, 39 D. J. Payne, 39 K. C. Schofield, 39 C. Touramanis, 39 A. J. Bevan, 40 
K. A. George, 40 F. Di Lodovico, 40 W. Menges, 40 R. Sacco, 40 G. Cowan, 41 H. U. Flaecher, 41 D. A. Hopkins, 41 
P. S. Jackson, 41 T. R. McMahon, 41 F. Salvatore, 41 A. C. Wren, 41 D. N. Brown, 42 C. L. Davis, 42 J. Allison, 43 
N. R. Barlow, 43 R. J. Barlow, 43 Y. M. Chia, 43 C. L. Edgar, 43 G. D. Lafferty, 43 T. J. West, 43 J. I. Yi, 43 C. Chen, 44 
W. D. Hulsbergen, 44 A. Jawahery, 44 D. A. Roberts, 44 G. Simi, 44 G. Blaylock, 45 C. Dallapiccola, 45 S. S. Hertzbach, 45 
X. Li, 45 T. B. Moore, 45 E. Salvati, 45 S. Saremi, 45 R. Cowan, 46 G. Sciolla, 46 S. J. Sekula, 46 M. Spitznagel, 46 
F. Taylor, 46 R. K. Yamamoto, 46 H. Kim, 47 S. E. Mclachlin, 47 P. M. Patel, 47 S. H. Robertson, 47 A. Lazzaro, 48 
V. Lombardo, 48 F. Palombo, 48 J. M. Bauer, 49 L. Cremaldi, 49 V. Eschenburg 49 R. Godang, 49 R. Kroeger, 49 
D. A. Sanders, 49 D. J. Summers, 49 H. W. Zhao, 49 S. Brunet, 50 D. Cote, 50 M. Simard, 50 P. Taras, 50 F. B. Viaud, 50 



2 



H. Nicholson, 51 N. Cavallo, 52 '0 G. De Nardo, 52 F. Fabozzi, 52 -© C. Gatto, 52 L. Lista, 52 D. Monorchio, 52 
P. Paolucci, 52 D. Piccolo, 52 C. Sciacca, 52 M. A. Baak, 53 G. Raven, 53 H. L. Snoek, 53 C. P. Jessop, 54 J. M. LoSecco, 54 
G. Benelli, 55 L. A. Corwin, 55 K. K. Gan, 55 K. Honscheid, 55 D. Hufnagel, 55 H. Kagan, 55 R. Kass, 55 J. P. Morris, 55 

A. M. Rahimi, 55 J. J. Regensburger, 55 R. Ter-Antonyan, 55 Q. K. Wong, 55 N. L. Blount, 56 J. Brau, 56 R. Frey, 56 
O. Igonkina, 56 J. A. Kolb, 56 M. Lu, 56 R. Rahmat, 56 N. B. Sinev, 56 D. Strom, 56 J. Strube, 56 E. Torrence, 56 A. Gaz, 57 

M. Margoni, 57 M. Morandin, 57 A. Pompili, 57 M. Posocco, 57 M. Rotondo, 57 F. Simonetto, 57 R. Stroili, 57 C. Voci, 57 
E. Ben-Haim, 58 H. Briand, 58 J. Chauveau, 58 P. David, 58 L. Del Buono, 58 Ch. de la Vaissiere, 58 O. Hamon, 58 

B. L. Hartfiel, 58 Ph. Leruste, 58 J. Malcles, 58 J. Ocariz, 58 A. Perez, 58 J. Prendki, 58 L. Gladney, 59 M. Biasini, 60 
R. Covarelli, 60 E. Manoni, 60 C. Angelini, 61 G. Batignani, 61 S. Bettarini, 61 G. Calderini, 61 M. Carpinelli, 61 

R. Cenci, 61 F. Forti, 61 M. A. Giorgi, 61 A. Lusiani, 61 G. Marchiori, 61 M. A. Mazur, 61 M. Morganti, 61 N. Neri, 61 
E. Paoloni, 61 G. Rizzo, 61 J. J. Walsh, 61 M. Haire, 62 J. Biesiada, 63 P. Elmer, 63 Y. P. Lau, 63 C Lu, 63 J. Olsen, 63 
A. J. S. Smith, 63 A. V. Telnov, 63 E. Baracchini, 64 F. Bellini, 64 G. Cavoto, 64 A. D'Orazio, 64 D. del Re, 64 E. Di 
Marco, 64 R. Faccini, 64 F. Ferrarotto, 64 F. Ferroni, 64 M. Gaspero, 64 P. D. Jackson, 64 L. Li Gioi, 64 M. A. Mazzoni, 64 
S. Morganti, 64 G. Piredda, 64 F. Polci, 64 F. Renga, 64 C Voena, 64 M. Ebert, 65 H. Schroder, 65 R. Waldi, 65 T. Adye, 66 
G. Castelli, 66 B. Franek, 66 E. O. Olaiya, 66 S. Ricciardi, 66 W. Roethel, 66 F. F. Wilson, 66 R. Aleksan, 67 S. Emery, 67 
M. Escalier, 67 A. Gaidot, 67 S. F. Ganzhur, 67 G. Hamel de Monchenault, 67 W. Kozanecki, 67 M. Legendre, 67 
G. Vasseur, 67 Ch. Yeche, 67 M. Zito, 67 X. R. Chen, 68 H. Liu, 68 W. Park, 68 M. V. Purohit, 68 J. R. Wilson, 68 
M. T. Allen, 69 D. Aston, 69 R. Bartoldus, 69 P. Bechtle, 69 N. Berger, 69 R. Claus, 69 J. P. Coleman, 69 M. R. Convery, 69 

J. C. Dingfelder, 69 J. Dorfan, 69 G. P. Dubois-Felsmann, 69 D. Dujmic, 69 W. Dunwoodie, 69 R. C. Field, 69 
T. Glanzman, 69 S. J. Gowdy, 69 M. T. Graham, 69 P. Grenier, 69 V. Halyo, 69 C. Hast, 69 T. Hryn'ova, 69 W. R. Innes, 69 
M. H. Kelsey, 69 P. Kim, 69 D. W. G. S. Leith, 69 S. Li, 69 S. Luitz, 69 V. Luth, 69 H. L. Lynch, 69 D. B. MacFarlane, 69 

H. Marsiske, 69 R. Messner, 69 D. R. Muller, 69 C. P. O'Grady, 69 V. E. Ozcan, 69 A. Perazzo, 69 M. Perl, 69 
T. Pulliam, 69 B. N. Ratcliff, 69 A. Roodman, 69 A. A. Salnikov, 69 R. H. Schindler, 69 J. Schwicning, 69 A. Snyder, 69 
J. Stclzer, 69 D. Su, 69 M. K. Sullivan, 69 K. Suzuki, 69 S. K. Swain, 69 J. M. Thompson, 69 J. Va'vra, 69 N. van Bakel, 69 
A. P. Wagner, 69 M. Weaver, 69 W. J. Wisniewski, 69 M. Wittgen, 69 D. H. Wright, 69 H. W. Wulsin, 69 A. K. Yarritu, 69 
K. Yi, 69 C. C. Young, 69 P. R. Burchat, 70 A. J. Edwards, 70 S. A. Majewski, 70 B. A. Petersen, 70 L. Wilden, 70 
S. Ahmed, 71 M. S. Alam, 71 R. Bula, 71 J. A. Ernst, 71 V. Jain, 71 B. Pan, 71 M. A. Saeed, 71 F. R. Wappler, 71 
S. B. Zain, 71 W. Bugg, 72 M. Krishnamurthy, 72 S. M. Spanier, 72 R. Eckmann, 73 J. L. Ritchie, 73 C. J. Schilling, 73 

R. F. Schwitters, 73 J. M. Izen, 74 X. C. Lou, 74 S. Ye, 74 F. Bianchi, 75 F. Gallo, 75 D. Gamba, 75 M. Pelliccioni, 75 
M. Bomben, 76 L. Bosisio, 76 C. Cartaro, 76 F. Cossutti, 76 G. Delia Ricca, 76 L. Lanceri, 76 L. Vitale, 76 V. Azzolini, 77 
N. Lopez-March, 77 F. Martinez- Vidal, 77 A. Oyanguren, 77 J. Albert, 78 Sw. Banerjee, 78 B. Bhuyan, 78 K. Hamano, 78 
R. Kowalewski, 78 I. M. Nugent, 78 J. M. Roney, 78 R. J. Sobie, 78 J. J. Back, 79 P. F. Harrison, 79 T. E. Latham, 79 
G. B. Mohanty, 79 M. Pappagallo, 79 ! H. R. Band, 80 X. Chen, 80 S. Dasu, 80 K. T. Flood, 80 J. J. Hollar, 80 
P. E. Kutter, 80 B. Mellado, 80 Y. Pan, 80 M. Pierini, 80 R. Prepost, 80 S. L. Wu, 80 Z. Yu, 80 and H. Neal 81 

(The BABAR Collaboration) 

1 Laboratoire de Physique des Particules, IN2P3/CNRS et Universite de Savoie, F-74941 Annecy-Le-Vieux, France 
Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 
3 Universitd di Bari, Dipartimento di Fisica and INFN, 1-70126 Bari, Italy 
* Institute of High Energy Physics, Beijing 100039, China 
5 University of Bergen, Institute of Physics, N-5007 Bergen, Norway 
6 Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 
8 Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44780 Bochum, Germany 
University of Bristol, Bristol BS8 1TL, United Kingdom 
1} University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 
11 Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 
12 Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 
13 University of California at Irvine, Irvine, California 92697, USA 
1Jt University of California at Los Angeles, Los Angeles, California 90024, USA 

15 University of California at Riverside, Riverside, California 92521, USA 

16 University of California at San Diego, La Jolla, California 92093, USA 

17 University of California at Santa Barbara, Santa Barbara, California 93106, USA 
18 University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 
19 California Institute of Technology, Pasadena, California 91125, USA 
20 University of Cincinnati, Cincinnati, Ohio 45221, USA 



3 



21 University of Colorado, Boulder, Colorado 80309, USA 
22 Colorado State University, Fort Collins, Colorado 80523, USA 
23 Universitat Dortmund, Institut fiir Physik, D-44221 Dortmund, Germany 
4 Technische Universitat Dresden, Institut fiir Kern- und Teilchenphysik, D-01062 Dresden, Germany 
25 Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, F-91128 Palaiseau, France 
26 University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 
27 Universita di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

28 Laboratori Nazionali di Frascati dell'INFN, 1-00044 Frascati, Italy 
29 Universita di Genova, Dipartimento di Fisica and INFN, 1-16146 Genova, Italy 
30 Harvard University, Cambridge, Massachusetts 02138, USA 
31 Universitat Heidelberg, Physikahsches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 
32 Imperial College London, London, SW7 2AZ, United Kingdom 
33 University of Iowa, Iowa City, Iowa 52242, USA 
34 Iowa State University, Ames, Iowa 50011-3160, USA 
' 5 Johns Hopkins University, Baltimore, Maryland 21218, USA 
36 Universitat Karlsruhe, Institut fiir Experimentelle Kernphysik, D-76021 Karlsruhe, Germany 
37 Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 ORSAY Cedex, France 
38 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
39 University of Liverpool, Liverpool L69 7ZE, United Kingdom 
40 Queen Mary, University of London, El 4NS, United Kingdom 
41 University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 
42 University of Louisville, Louisville, Kentucky 40292, USA 
43 University of Manchester, Manchester M13 9PL, United Kingdom 

44 University of Maryland, College Park, Maryland 20742, USA 
4 '' University of Massachusetts, Amherst, Massachusetts 01003, USA 
46 Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 
4? McGill University, Montreal, Quebec, Canada H3A 2T8 
48 Universita di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 
49 University of Mississippi, University, Mississippi 38677, USA 
50 Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3J7 
51 Mount Holyoke College, South Hadley, Massachusetts 01075, USA 
52 Universita di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 
53 NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

54 University of Notre Dame, Notre Dame, Indiana 46556, USA 
55 Ohio State University, Columbus, Ohio 43210, USA 
56 University of Oregon, Eugene, Oregon 97403, USA 
57 Universita di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 
58 Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 
Universite Denis Diderot- Paris7, F-75252 Paris, France 
59 University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 
60 Universita di Perugia, Dipartimento di Fisica and INFN, 1-06100 Perugia, Italy 
61 Universita di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56127 Pisa, Italy 
62 Prairie View A&M University, Prairie View, Texas 77446, USA 
Princeton University, Princeton, New Jersey 08544> USA 
4 Universita di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 
65 Universitat Rostock, D-18051 Rostock, Germany 
66 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 
67 DSM/Dapnia, CEA/Saclay, F- 9 1191 Gif-sur-Yvette, France 
68 University of South Carolina, Columbia, South Carolina 29208, USA 
6 ' Stanford Linear Accelerator Center, Stanford, California 94309, USA 
70 Stanford University, Stanford, California 94305-4060, USA 
71 State University of New York, Albany, New York 12222, USA 

72 University of Tennessee, Knoxville, Tennessee 37996, USA 

73 University of Texas at Austin, Austin, Texas 78712, USA 
74 University of Texas at Dallas, Richardson, Texas 75083, USA 

75 Universita di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 
76 Universita di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 
IFIC, Universitat de Valencia- C SIC, E-46071 Valencia, Spain 
78 University of Victoria, Victoria, British Columbia, Canada V8W 3P6 
79 Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 
80 University of Wisconsin, Madison, Wisconsin 53706, USA 
81 Yale University, New Haven, Connecticut 06511, USA 



4 



(Dated: February 3, 2008) 

We report the results of a CP violation analysis of the decay B ± — ► D 7V + jr - 7V oK ± , where D^+^-^o 
indicates a neutral D meson detected in the final state 7r + 7r _ 7r°, excluding K®n°. The analysis 
makes use of 324 million e + e~ — > BB events recorded by the BABAR experiment at the PEP-II 
e + e~ storage ring. By analyzing the ir + n~ir° Dalitz plot distribution and the B ± — + D^+^- ^oK^ 1 
branching fraction and decay rate asymmetry, we calculate parameters related to the phase 7 of 
the CKM unitarity triangle. We also measure the magnitudes and phases of the components of the 
D° — * 7T + 7r~7r° decay amplitude. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



An important component of the program to study 
CP violation is the measurement of the angle 7 = 
arg {—V ud V* b /V cd V* b ) of the unitarity triangle related 
to the Cabibbo-Kobayashi-Maskawa quark mixing ma- 
trix The decays B -> £>(*)° can be used to mea- 
sure 7 with essentially no hadronic uncertainties, exploit- 
ing interference between b — > ucs and b — ► cus decay am- 
plitudes [2] . In one of the measurement methods 0] , 7 is 
extracted by analyzing the £>-decay Dalitz plot distribu- 
tion in B ± — > DK ± with multi-body D decays 0. This 
method has only been used with the Cabibbo-favored de- 
cay D — > K®it +/ k~ @,@|, and Cabibbo-suppressed decays 
are expected to be similarly sensitive to 7 [7| . We present 
here the first CP-violation study of B ± — > DK with a 
multibody, Cabibbo-suppressed D decay, D — > h + it~it . 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II e + e~ storage ring, running 
on the T(4S) resonance. Samples of simulated Monte 
Carlo (MC) events were analyzed with the same recon- 
struction and analysis procedures. These samples include 
an e + e~ — > BB sample about five times larger than the 
data; a continuum e + e~ — > qq sample, where q is a u, 
d, s, or c quark, with luminosity equivalent to the data; 
and a signal sample about 300 times larger than the data, 
with both phase space D decays and decays generated ac- 
cording to the amplitudes measured by CLEO §|. The 
BABAR detector and the methods used for particle recon- 
struction and identification are described in Ref. 

The reader is referred to Ref. [13] for details of the 
event selection criteria. Briefly, we use event-shape vari- 
ables to suppress the continuum background, and iden- 
tify kaon and pion candidates using specific ionization 
and Cherenkov radiation. The invariant mass of D can- 
didates must satisfy 1830 < M D < 1895 MeV/c 2 . We 
require 5272 < ities < 5300 MeV/c 2 , where tues = 
\/E^ M /4: — |pb| 2 , Eqm is the total e + e~ center-of-mass 
(CM) energy, and pb is the B candidate CM momen- 
tum. Events must satisfy -70 < AE < 60 MeV, 
where AE = Eb — Eqm/ '2 and Eg is the B candidate 
CM energy. We exclude the decay mode D — » A^Tr , 
which is a previously studied CP eigenstate not related 
to the method of Ref. Q, by rejecting candidates with 
489 < M(7r+7T~) < 508 MeV/c 2 or for which the dis- 
tance between the it + tt~ vertex and the B~ candidate 
decay vertex is more than 1.5 cm. We reject B^ — > 



D^+^-^oK candidates in which the K if^ invariant 
mass satisfies 1840 < M(K ± n^ z ) < 1890 MeV/c 2 , to 
suppress B~ — > -D^-^+P" decays. We require d > 0.25, 
where d [13] is a neural net variable that separates signal 
candidates (which peak toward d = 1) from those with 
a misreconstructed D (peaking toward d = 0). In events 
with multiple candidates, we keep the candidate whose 
tties value is closest to the nominal B^ mass [Tl| . 

For each B ± — ► D 7I +^- v oK ± candidate, we compute 
the neural net variable q [lflj • The q distribution of BB 
events peaks toward q = 1, while that of continuum 
peaks at q = 0. For v S {q,d}, we define the variables 
v' = tanh -1 [(1/ - \{v max + v mln ))l\(y max + v min )] , 
where q max = d max = 1, q min = 0.1, and d mm = 0.25 
are the allowed ranges for q and d. The v' variables can 
be conveniently fit with Gaussians, as described later. 

As in Ref. [l(| , we identify in the MC samples ten event 
types, one signal and nine different backgrounds. We list 
them here with the labels used to refer to them through- 
out the paper. DK s i s : B^ — » D 7T +^-^oK ± events that 
are correctly reconstructed; these are the only events 
considered to be signal. DK hsd : B ± — » D^+^-^o A± 
events that are misreconstructed; namely, some of the 
particles used to form the final state do not originate 
from the B ± — ► D v + n - w oK decay. D-kd {Dirp): 
B~ — > D°ir~, D° — > 7T + 7r~7r decays, where the de- 
cay D° — > 7T + 7r _ 7r is correctly reconstructed (misrecon- 
structed). DKX. B -> £>(*)#(*)- events not contain- 
ing the decay D — ► 7r + 7r~7r°. DtvX: B — » D^n~ and 
fl — > D^p~ decays, excluding 13 — » 7r + 7r~7r°. BBCd 
(BBCp): all other BB events with a correctly recon- 
structed (misreconstructed) D candidate, qqo (qqp): 
continuum e + e~ — ► qq events with a correctly recon- 
structed (misreconstructed) D candidate. 

The measurement of the CP parameters proceeds in 
three steps, each involving an unbinned maximum likeli- 
hood fit. In step 1, we measure the complex Dalitz plot 
amplitude ct(s+, s_) for the decay D° — > 7r + 7r~7r°, where 
s± — m 2 (ir 7T°) are the squared invariant masses of the 
tt^it pairs. In step 2, we extract the numbers of B + and 
B~ signal events and background yields. We obtain the 
CP parameters in step 3. 

We parameterize a(s+,s_) using the isobar model, 
«(«+,*-) = [aNRe #NR + £ r a r e l ^A r (s + ,s_)]/A a , 
where the first term represents a nonresonant contri- 
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bution, the sum is over all intermediate two-body reso- 
nances r, and N a is such that J ds + ds-\a(s + , s-)| 2 = 1. 
The amplitude for the decay chain D° — > rC, r — > Ai? 
is A r (s+, s_) = F r F s (m 2 . - - m r r r (M AB )) _1 , 
where m r is the peak mass of the resonance M\ B 
is the squared invariant mass of the AB pair, F r is a 
spin-dependent form factor 12], and r r (M4s) is the 
mass-dependent width for the resonance r |12j | . The spin 
factors F s are i*o = m 2 D , F 1 — M BC — M 2 AC + {m 2 D — 

m2 c)( m A ~ ™b) m ab> and F 2 = (if -\^ cd ^ab) m D 2 ' 
where ^ = M\ B - 2m 2 - 2ml + (m 2 - ml) M~ k 2 , and 
rrii is the mass of particle i . 

In step 1, we determine the parameters onr, a r , 
</>nr, and 4> r by fitting a large sample of -D and D° 
mesons, flavor-tag ged through their production in the de- 
cay D* + — > D°7r + [13j]. To select this sample, we require 
the CM momentum of the D* candidate to be greater 
than 2770 MeV/c, and \M D , - M D - 145.4 MeV/c 2 | < 
0.6 MeV/c 2 , where Md* is the invariant mass of the D* 
candidate. The signal and background yields are ob- 
tained from a fit to the M B distribution, modeling the 
signal as a Gaussian and the background as an exponen- 
tial. The signal Gaussian peaks at 1863.7 ± 0.4 MeV/c 2 
and has a width of 17.4 ± 0.8 MeV/c 2 . 

Of the D° candidates in the signal region 1848 < 
M D < 1880 MeV/c 2 , we obtain from the fit N s = 
44780 ±250 signal and N B = 830 ±70 background events. 
To obtain the parameters of a(s±, s^), we fit these can- 
didates with the probability distribution function (PDF) 
N s |a(s + ,s_)| 2 e(s + ,s_) + N B \f B (s+, s-)| 2 , where the 
background PDF /s(s_|_,s_) is a binned distribution 
obtained from events in the sideband 1930 < Md < 
1990 MeV/c 2 , and e(s+,s_) is an efficiency function, 
parameterized as a two-dimensional third-order polyno- 
mial determined from MC. To within the MC-signal sta- 
tistical uncertainty, e(s+,s_) = e(s_,s+). The region 
M D < 1848 MeV/c 2 , which contains D° -> R-tt+tt 
events that are absent from the signal region, is not used. 

Table U summarizes the results of this fit, with system- 
atic errors obtained by varying the masses and widths of 
the p(1700) and a resonances, setting F r = 1, and vary- 
ing e(s_|_,s_) to account for uncertainties in reconstruc- 
tion and particle identification. The Dalitz plot distribu- 
tion of the data is shown in Fig. [TJa-c). The distribution 
is marked by three destructively interfering pn ampli- 
tudes, suggesting an / = 0-dominated final state [T3 |. 

The fit for step i G {2, 3} uses the PDF 

where & is the set of Hi event variables £i = {AE, q' , d'}, 
£,2 = {AE, q' , s_, s + }, t corresponds to one of the ten 
event types listed above, N t = + N[~ is the number 
of events of type t, A t = (Nf — N*)/N t is their charge 
asymmetry, C = ±1 is the electric charge of the B can- 



TABLE I: Result of the fit to the D* + -> D°tt + sam- 
ple, showing the amplitudes ratios R r = o, r /a p +(rro)> phase 
differences A<j> r = <f> r — <f> P +(no)i an d fit fractions f r = 
J \a r A r (s+, s~)\ 2 ds~ds+. The first (second) errors are statis- 
tical (systematic). We take the mass (width) of the a meson 
to be 400 (600) MeV/c 2 . 



State 


Rr (%) 


A<p r (") 


M%) 


p+(770) 


100 





67.8±0.0±0.6 


p°(770) 


58.8±0.6±0.2 


16.2±0.6±0.4 


26.2±0.5±1.1 


P~(770) 


71.4±0.8±0.3 


-2.0±0.6±0.6 


34.6±0.8±0.3 


p+(1450) 


21±6±13 


-146±18±24 


0.11±0.07±0.12 


p°(1450) 


33±6±4 


10±8±13 


0.30±0.11±0.07 


p-(1450) 


82±5±4 


16±3±3 


1.79±0.22±0.12 


p+(1700) 


225±18±14 


-17±2±3 


4.1±0.7±0.7 


p°(1700) 


251±15±13 


-17±2±2 


5.0±0.6±1.0 


p~(1700) 


200±11±7 


-50±3±3 


3.2±0.4±0.6 


/o(980) 


1.50±0.12±0.17 


-59±5±4 


0.25±0.04±0.04 


/o(1370) 


6.3±0.9±0.9 


156±9±6 


0.37±0.11±0.09 


/o(1500) 


5.8±0.6±0.6 


12±9±4 


0.39±0.08±0.07 


/o(1710) 


11.2±1.4±1.7 


51±8±7 


0.31±0.07±0.08 


/ 2 (1270) 


104±3±21 


-171±3±4 


1.32±0.08±0.10 


cr(400) 


6.9±0.6±1.2 


8±4±8 


0.82±0.10±0.10 


Non-Res 


57±7±8 


-11±4±2 


0.84±0.21±0.12 



didate, and rj = ^2 t N t . Using MC, we verify that the 

£i and £j (i ^ j) distributions are uncorrelated for each 

(c) 

event type. Therefore, the PDFs V\ t arc the products 

V 2 , t (AE,q',d') = £ t (AE)Q t (q , )C t (d') 
VC t (AE,q', s+ ,s-) = £ t {AE)Q t {q')V' C t { s+ ,s^).{2) 

The parameters of the Dalitz plot PDF T> DK Bis ( s +- s -) 
are obtained from the data as described below. Those 
of all other functions in Eq. are obtained from the 
MC samples. The functions £t(AE) are parameterized 
as the sum of a Gaussian and a second-order polyno- 
mial. The PDFs Qt(q') and C t (d') are the sum of a 
Gaussian and an asymmetric Gaussian. The PDF pa- 
rameters are different for each event type. Assuming no 
CP violation in the background, we take V'f (s + , S-) = 
X>7(s_,s+) and A t = for t ^ DK sig . The func- 
tions V M (s|, s_) and D' DKh d (s_|_, s_) are binned his- 
tograms obtained from the MC. For other event types, 
V t (s + , s_) = e(s + , s_)Vf (s + , s_), where the efficiency 
function e(s + ,s_) has different parameters for well- 
reconstructed and misreconstructed D candidates. 

The signal Dalitz PDF accounts for interference be- 
tween the b — > ucs and b — > cus amplitudes A u and A c : 

P M, iE ( s +' s -) = K S T' S ±) + z±a(s±,s T )\ 2 , (3) 

where z± = |A u /^4 c |e^' 5±7 ^ and S is a CP-even phase. 

In the step-2 fit, we extract the B ± — » D^+„- n o 
signal yield and asymmetry, as well as some background 
yields, as described in Ref. 10]. From this fit we find 
NDK Bi — 170 ± 29 signal events and a decay rate asym- 
metry Ar>K Big = —0.02 ±0.15. Errors are statistical only. 
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Only the complex parameters z± are free in the step-3 
fit. This fit minimizes the function 



^ 1 



e=l 



(4) 



where iV ev is the number of events in the data sample. 
The term \ 2 — J2u v =i ^uV~^X v increases the sensitiv- 
ity of the fit by using the results of the step-2 fit via 



where 



x 2 



n± 



N D k As -(n-+n+), 

A D k sis - (n- - + n+), (5) 



/ |a(s T , s±)| 2 e(s + , s-)ds+ds- 



(6) 



are the expected numbers of -E^ signal events. In 
Eq. (JBJ), A^ is the product of the number N B + B - of 
charged B + B~ pairs in the dataset, the branching frac- 
tions B(B~ -► D°K~) [IH and B{D° -> tt+tt"^ ) [Hj], 
and the total reconstruction efficiency e = 11.4%. The 
error matrix V uv is the sum of two components: the 
step-2 fit error matrix T^* at , which is almost diagonal 
(the correlation coefficient is —2.8%), and the ./V sys- 
tematic error matrix V^ st . Here V^" 5 * = = 0, and 
V^l — J2c=i(N° c^ 1 ) 2 ; where cr£ cl are the relative er- 
rors on the four components N B + B - (1.1%), e (3.3%), 
B(D -> tt+tt-tt ) (3.8%) H3, and B(B~ -» D°K~) 



(5.9%) [ll|. 

We parameterize z± with the polar coordinates 

P± ,\ Z± - X0 \, fl±ste n-i(_^_), (7) 

where xo is a coordinate transformation parameter, 

x o = -J & M>+, s_) a* (s_, s+)] ds + ds_ = 0.850. (8) 

This parameterization is optimal due to the polar sym- 
metry of n± = N°(l + p± 2 — xq 2 ). Other parameteri- 
zations, such as (\A u /A c \,j,S) or (3?[z±], 3?[z±]), result 
in significant nonlinear correlations between the fit vari- 
ables, which cannot be parameterized with an error ma- 
trix, and bias the fit result. The polar coordinates enable 
a significant improvement in sensitivity due to the x 2 
term in Eq. ([4J, and are determined from parameterized 
simulation to be unbiased. The step-3 fit yields 

p- = 0.72 ± 0.11 ± 0.04, 6L = (173 ± 42 ± 2)°, 

P+ = 0.75 ±0.11 ±0.04, 6»+ = (147 ±23 ± 1)°, (9) 

where the first errors are statistical and the second are 
systematic, due only to Vlf st . The largest correlation 
coefficient is c p _ p+ = 14%, originating from V^f st . All 
others are 1% or less. Contours of constant C values are 
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FIG. 1: (Charge conjugation is implied for all plots.) (a,b) 
Projections of the D* + — > D°tv + data events and PDF onto 
the Dalitz plot variables s+ and s_. (c) The 2- dimensional 
(s+,s_) distribution of the D* + — * D°tt + data, (d) One-, 
two-, and three-standard-deviation contours of C as a func- 
tion of 9± vs. p±. The solid (dashed) curves correspond to 
B + (B~) results. The no-interference point (p± = xo,0± = 
180°) is marked with an x. (e,f) Projection of the B~ — » 
D^+^-^oK" candidate data onto s + and s_. 



shown in Fig. QJd). Projections of the data and the PDF 
onto s+ and s_ are shown in Fig. [ljc-f). 

Additional systematic errors due to the analysis pro- 
cedure are evaluated for the signal branching fraction, 
charge asymmetry, p±, and 9±. The uncertainty in the 
model used for a(s + ,s_) is the largest source of error 
on the CP parameters: cr p n ± odcl = 0.03, af° del = 14°, 
^modoi — i\o rp n j s error i s evaluated by removing all but 
the p(770), p(1450), /o(980), and nonresonant terms in 
a(s+, s_); adding an /2(1525), an to, and a nonresonant 
P-wave contribution; varying the meson "radius" param- 
eter in F r [HI; and propagating the errors from Tabled] 
Uncertainties due to the masses and widths of the p(1700) 
and a resonances are small by comparison. Other errors 
are due to uncertainties on background yields that are 
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fixed in the fits [lOJ], finite MC sample size, a possible re- 
construction efficiency charge asymmetry, and uncertain- 
ties in the background PDF shapes, evaluated by com- 
paring MC and data in signal-free sidebands of the vari- 
ables Md, AE, and itles- We also evaluate errors due to 
possible charge asymmetries in DKX and -Diibgd, uncer- 
tainties in particle identification and the efficiency func- 
tions, the finite s± measurement resolution, the back- 
ground PDF Jb in the D* sample, D-flavor mistagging 
in the D* sample, and correlations between the D flavor 
and the kaon charge in qqu events. These errors add in 
quadrature to a s /f = 0.05, a s g yst = 19°, er^ st = 13°, and 
are combined with the systematic errors of Eqs. (J9j> . 

The analysis procedure is validated in several ways. 
Conducting the analysis on the MC sample yields results 
consistent with the generated values. We carry out the 
step-3 fit on a sample of 1800 ± 70 B~ -> I^+^oTr" 
events, obtaining the background Dalitz plot distribution 
from the AE sideband. The fit yields p- = 0.815 ±0.034, 
9- = (186±7)°, p+ = 0.854±0.035, 9+ = (192±7)°, con- 
sistent with p± = xq, 9± = 180°, which corresponds to 
z± = 0. We verify the signal efficiency by measuring the 
branching fraction B(B~ -> D°tt-) with D° -> R-tt+tt 
and D° — > 7r + 7r~7r°. We compare the fit variable distri- 
butions of data and MC events in signal- free sidebands. 
Good agreement is found in all cases. 

In summary, using a sample of (324.0 ± 3.6) x 10 6 
e + e _ — > BB events, we observe 170 ± 29 B — > 
D^+^-^aK^ events. We calculate the branching fraction 
and decay rate asymmetry 

B{B ± -> D^+k-^K^ = (4.6 ±0.8 ±0.7) x 10~ 6 , 
A(B ± ->• D^-^K^) = -0.02 ±0.15 ±0.03, (10) 

and the CP- violation parameters 

p- = 0.72 ± 0.11 ± 0.06, 9- = (173 ± 42 ± 19)°, 

p+ = 0.75 ± 0.11 ± 0.06, 6+ = (147 ± 23 ± 13)°, (11) 

where the first errors are statistical and the second are 
systematic. The parameters p±, 9± are defined in Eq. 10. 
While the errors on 9± are too large for a meaningful de- 
termination of 7 with these results alone, our errors on 
p± are small enough to make a non-negligible contribu- 
tion to the overall precision of 7 in a combination of all 
measurements related to 7. In addition, we measure the 
magnitudes and phases of the components of the ampli- 
tude of the decay D° — » Tr + ir~ir° in the isobar model. 
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